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REVERSE  CURRENT  IN  SOLAR  FLARES 
by 

J.W.  Knight  and  P.A.  Sturrock 

Institute  for  Plasma  Research 
and 

Department  of  Applied  Physics 
Stanford  University 
Stanford,  California 

ABSTRACT 

We  examine  the  proposal  that  impulsive  x-ray  bursts  are  produced 
by  high-energy  electrons  streaming  from  the  corona  to  the  chromosphere. 

We  find  that  the  currents  associated  with  these  streams  are  so  high 
that  either  the  streams  do  not  exist  or  their  current  is  neutralized 
by  a reverse  current. 

Analysis  of  a simple  model  indicates  that  the  primary  electron 
stream  leads  to  the  development  of  an  electric  field  in  the  ambient 
corona  which  (a)  decelerates  the  primary  beam  and  (b)  produces  a 
neutralizing  reverse  current.  It  appears  that,  in  some  circumstances, 
this  electric  field  could  prevent  the  primary  beam  from  reaching  the 
chromosphere.  In  any  case,  the  electric  field  acts  as  an  energy  exchange 
mechanism,  extracting  kinetic  energy  from  the  primary  beam  and  using  it 
to  heat  the  ambient  plasma.  This  heating  is  typically  so  rapid  that  it 
must  be  expected  to  have  important  dynamical  consequences. 
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I . Introduction 


Much  of  the  line  emission  from  solar  flares  must  originate  in 
plasma  of  chromospheric  density.  On  the  other  hand,  the  corona  offers 
a much  larger  volume  for  the  storage  of  magnetic  energy  than  does  the 
chromosphere.  .’his  has  led  to  models  in  which  free  magnetic  energy, 
released  in  the  corona,  propagates  to  the  chromosphere  in  the  form  of 
streams  of  high-energy  particles  (Sweet,  1969;  Sturrock,  197^)*  Obser- 
vational evidence  for  high-energy  electron  streams  comes  from  analysis 
of  microwave  bursts  and  impulsive  x-ray  bursts  (Kane,  1972).  Evidence 
concerning  the  origin  of  impulsive  x-ray  bursts  is  conflicting  and  it 
is  possible  that  some  radiation  originates  at  the  chromosphere,  as 
analyzed  by  Brown  (1972),  Hudson  (1972)  and  Petrosian  (1973),  and  that 
some  originates  in  the  corona  as  proposed  by  de  Jager  and  Kundu  (1963), 
Kane  and  Anderson  (1970)  and  Kane  (1972). 

It  has  recently  occurred  to  us,  and  to  other  solar  physicists 
including  Brown  (1976)  and  Colgate  (1976),  that  the  electrical  currents 
associated  with  the  proposed  electron  streams  are  so  high  that  the 
streams  either  do  not  exist  or  are  neutralized  by  reverse  currents  in 
the  ambient  plasma.  The  reason  for  this  statement  are  discussed 
briefly  in  Section  II.  An  idealized  model  in  which  the  primary  current 
is  neutralized  by  a reverse  current  is  presented  in  Section  III  and 
applied  to  solar  flares  in  Section  IV.  Some  of  the  consequences  of 


these  considerations  are  discussed  in  Section  V. 


1 1 . Magne tic  Field  Associated  with  Proposed  Electron  Stream 


IE  energy  released  in  the  corona  propagates  to  the  chromosphere  in 
the  form  of  high-energy  electrons,  this  electron  stream  will  constitute 
an  electrical  current.  If  N is  the  total  number  of  electrons  streaming 
downward  over  the  duration  t (sec.)  of  the  impulsive  phase,  the 
magnitude  J (e.m.u.)  of  the  current  may  be  estimated  from 


J ~ e c'1  N T*1 


(2.1) 


If  the  transverse  and  longitudinal  dimensions  of  the  stream  are 
of  order  L (cm.),  an  estimate  of  the  strength  B (gauss)  of  the  magnetic 
field  produced  by  the  stream  is  given  by 


B ~ 2J  L 


and  the  total  energy  U (erg)  of  this  magnetic  field  may  be  estimated 


i_  T 3 b2  i_  2 
Btt  L B ~ 2n  J 


(2.3) 


which  becomes 


2 -2  m2  -2  T 

e c N T L 


-4o.4  .,2  -2 


(2.4  ) 


Kane  and  Anderson  (1970)  estimate  the  total  energy  involved  in  a 

29  2 

typical  small  flare  to  be  ~ 10  erg,  the  time  scale  to  be  ~ 10  sec, 

8.5 

the  characteristic  length  scale  to  be  ~ 10  * cm,  and  infer  from  the 

35 

x-ray  data  that  the  total  number  of  energetic  electrons  is  ~ 10  . For 

13-2  5 

these  values,  the  above  formulas  lead  to  estimates  J ~ 10  , B ~ 10 

and  U ~ 10^  . It  is  clear  that  this  model,  involving  energy  transport  by 
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an  unneutralized  electron  beam,  leads  to  unacceptably  high  values  of  the 
magnetic  field  of  the  beam  and  magnetic  energy  associated  wit  1 the  beam. 
When  one  considers  large  flares  rather  than  small  flares,  one  faces  the 
same  difficulty. 

One  resolution  of  this  difficulty  is  to  abandon  the  assumption  that 
energy  is  transported  from  the  corona  to  the  chromosphere  by  an  electron 
stream.  One  might  assume,  for  instance,  that  energy  is  transported  as 
heat  flux.  The  above  parameters  for  a small  flare  then  indicate  that 
the  mean  heat  flux  F (erg  cm  ‘ s ^ ) would  be  ~ 10^  . From  the  formula 
appropriate  for  a fully  ionized  plasma  (Spitzer,  I962), 


F . „ f ~ 1C’6  T5/2  f , 
dz  dz 


(2.5) 


we  find  that  the  electron  temperature  T (°K)  must  be  at  least  10  to 

give  the  required  heat  flux.  For  the  concept  of  heat  flux  to  be  valid, 

the  mean  free  path  must  be  small  compared  with  the  length  scale  for 

-1  7 5 

temperature  variation,  e.g.  10  L,  which  is  ~ 10  * for  the  case  now 
being  considered.  Since  the  mean  free  path  y for  electron-electron 
and  electron-ion  collisions  is  given  by 

y « 10U  T2  n"1  , (2.6! 

we  see  that  the  electron  density  n (cm  must  be  at  least  lO^'^. 

This  is  substantially  higher  than  is  typical  in  the  normal  corona  at 
the  initiation  of  a flare,  although  densities  of  this  magnitude  do  build 
up  in  the  course  of  a flare.  This  indicates  that,  although  heat  con- 
duction may  be  an  important,  and  perhaps  the  dominant,  mechanism  of 


energy  transfer  in  the  later  stages  of  a flare,  it  cannot  be  the  principal 


mode  of  energy  transfer  at  flare  onset.  (One  should,  however,  bear  in 
mind  the  possibility  that  instabilities  may  lead  to  modification  of  the 
transport  properties  of  the  plasma,  modifying  both  the  heat  conduction 
and  effective  mean  free  path.) 
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I * 


III.  Model  of  Electron  Stream  and  Reverse  Current: 

It  is  well  known  that  a plasma  tends  to  preserve  charge  neutrality. 

A process  which  tends  to  give  an  excess  positive  or  negative  charge  in 
some  region  will  lead  to  electric  fields  acting  upon  the  plasma.  Movement 
of  electrons  in  response  to  this  electric  field  will  then  restore  charge 
neutrality.  One  expects  that  analogous  processes  will  also  tend  to 
maintain  current  neutrality.  If  an  electron  beam  is  suddenly  introduced 
into  a olasma,  there  is  a sudden  change  in  the  magnetic  field  structure 
which  will  develop  Induced  electric  fields  opposing  the  primary  current. 
For  this  reason,  we  should  investigate  the  possibility  that  a reverse 
current  develops  in  the  corona  balancing  (or  almost  balancing)  the 
primary  current  due  to  an  assumed  electron  stream  flowing  from  the 
corona  to  the  chromosphere. 

Our  considerations  are  substantially  simplified  if  we  consider 
a symmetrical  flux  tube  rooted  in  the  photosphere  and  extending  up 
into  the  corona,  and  assume  that  electrons  are  accelerated  at  the  top 
of  the  flux  tube  by  the  development  of  stochastic  electric  fields  or 
by  some  other  mechanism.  The  ejection  of  these  electrons  down  towards 
the  chromosphere  then  leads  to  a charge  imbalance  which  leads  to  the 
development  of  an  electric  field  of  such  a magnitude  that  charge 
balance  will  be  restored.  If  there  is  no  net  change  in  charge  density, 
whatever  current  flows  along  tiie  flux  tube  must  be  constant  along  the 
length  of  the  tube.  If  the  tube  is  symmetrical,  this  current  must  be 
zero.  Hence  the  strong  tendency  of  plasma  to  remain  charge  neutral 
implies  that  any  primary  current  will  generate  a neutralizing  secondary 
reverse  current. 
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In  order  to  simplify  this  exploratory  analysis,  we  consider  only 
a steady-state  situation.  We  also  assume  that  the  background  plasma 
can  be  adequately  described  by  a Maxwellian  velocity  distribution  and 
use  transport  coefficients  based  on  this  assumption  (Spitzer,  1962  1 . 

We  are  interested  in  the  case  that  the  primary  electron  stream  is 
composed  of  high-energy  electrons  which  consequently  have  long  mean 
free  paths.  However,  we  shall  find  that  an  electric  field  davelops 
which  decelerates  the  electron  stream.  But  when  the  electron  energy 
becomes  comparable  with  the  thermal  energy,  r mean  free  path  will  be 
sufficiently  short  that  the  primary  electrons  will  merge  with  the  back- 
ground plasma.  As  a simple  representation  of  this  process,  we  ignore 
collisions  in  discussing  the  primary  beam  but  we  assume  that  an  electron 
of  the  primary  beam  is  absorbed  into  the  background  plasma  when  it  is 
decelerated  to  zero  energy. 

If,  as  a further  simplification,  we  consider  a magnetic  flux  tube 
of  constant  cross  section,  we  may  use  the  following  simple  one-dimensional 
form  of  the  Vlasov  equation: 


vdf  e dtg  9f 

9s  m ds  9v 


(3.1) 


where  s measures  arc  length  along  the  tube,  v is  velocity  (along  the  tube', 
f(s,v)  is  the  velocity  distribution  function  of  the  primary  electron 
stream,  and  cp  is  the  electrostatic  potential. 

At  the  top  of  the  loop  (s  = 0),  the  primary  electron  stream  is 
moving  with  positive  velocity  on  the  side  corresponding  to  increasing 
s.  Hence  we  may,  without  ambiguity,  express  f in  terms  of  defined  by 


2e 


w 


(3.2) 


The  initial  distribution  function  may  therefore  be  expressed  as 


f(o,  v)  = f(i|i)  . '3.3 

With  this  initial  condition,  we  find  that  the  solution  of  the  Vlasov 
equation  (3.1)  is 

f(s,  v)  = F (ill  - cp)  . (3.4) 


The  current  density  j (e.m.u.)  in  the  primary  electron  stream 


is  given  by 


= - f / f(s, 


")  vdv 


(3.5) 


which  may  be  expressed  as 


me 


/%(,- 

J0 


cp)  di|f  . 


(3.6) 


Since  cp  will  prove  to  be  negative  in  the  region  of  interest,  it 
is  convenient  to  write 


9 = - cp 


(3.  ') 


so  that  (3-6)  may  be  reexpressed  as 

2 


_ __a!  f 

s me  J . 


F(X)  dX 


(3.8) 


Charged  neutrality  requires  that 


j + j =0 
Jp  Js 


(3.9' 


where  j is  the  secondary  current  induced  in  the  background  plasma.  We 
here  assume  that  the  density  and  temperature  are  such  that  this  may  be 


i 
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represented  by  Ohm's  law. 


j = aE  = a t-  , 
Jp  ds 


where  <y  is  the  electrical  conductivity  expressed  in  modified  gaussian 
units . 

It  is  convenient  to  introduce  a new  independent  vairable  § to 
replace  s,  by  the  relationship 


df  = — r^T  • 

Then,  on  substituting  equations  (3*8)  and  (3. 10)  into  (3-9)  and 
differentiating  with  respect  to  f,  we  obtain 

^ + £F<9>3i-°  • 

d5?  me  d? 

It  is  convenient  to  solve  this  equation  for  5 in  terms  of  0, 

I = x(e)  , 

rather  than  vice  versa.  Equation  (3.13)  becomes 


me  /dX  _ . . 

jM  ^-r(e)  • 


which  may  be  integrated  once  to  give 


(3.11) 


(3.12 


(3.13) 


(3.15 


if  we  assume  that  9 = 0 (tp  = 0)  and  X = 0 ( ^ = 0 ) at  s = 0.  We  find, 
from  (3.8),  (3.9)  and  (3.I0),  that 
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IV.  Specific  Model 


It  is  now  convenient  to  introduce  a specific  form  for  F > \|i  ) , 


F (if)  = K(iji0  + dr)’ 


-.1 


This  is  a power-law  distribution  at  high  energy  which  flattens  at  low 
energy,  the  "knee"  being  characterized  by  i|/  . 

*2  - 1 

We  introduce  the  symbol  H(i|f,  s)  for  the  flux  of  electrons  (cm  s ' 


of  energy  exceeding  et|f  at  the  position  s 


y.oo 

F(*  + 0(b))  df  • 

lit 


(4.2' 


If  the  initial  flux  is  written  as  H (i(f),  we  find  that 


V*>-  <*o  + *>'Y  + 1 


(4.3) 


so  that  the  total  particle  flux  is  given  by 


V0)  = FihrC+1 


(4.4) 


With  the  form  (4.1)  for  F ( \|» ) , equation  (3.17)  integrates  to  give 


x(e ) 


Y-l  me 
V e?K 


(V  + 9)v  - 


(4.5) 


We  easily  obtain  from  (4.5)  an  expression  for  the  (negative)  electric 
potential  9 in  terms  of  the  modified  distance  measure  §: 


e(S)  = 


+ efK  5 

^0  y-l  me 


1/Y 


*0  * 


(4.6) 


Hence,  from  (4.2),  we  find  that 


h(v.  s) 


eK 


(Y“ 1 ) m 


Ly  + _2£-  fifK  .1 


1/Y‘ 


til  + Ml ' + — L 

y > Tn  v- 


*0  T Y-l  me 


(4.7) 
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t 


On  noting  that  the  electric  current  carried  by  the  stream  is 
related  to  H(t)i,  s)  by 


j „(s)  = - 7 H (0 , s)  , 


we  see  that 


Vs  j 


e2K  y _ 

(y- 1 ) me  x>  y 


2 "I 

_5L.  O. 

Y-l  me  s 


In  order  to  specify  the  current,  particle  flux  and  electric  field 
as  functions  of  s,  we  must  adopt  a specific  form  for  a(s).  A convenient 
approximation  to  the  structure  of  the  solar  atmosphere,  which  is 
expressible  in  analytic  form,  is  provided  by  the  constant  heat  flux 
model.  If  we  now  assume  that  s measures  distance  vertically  downwards 
from  the  coror.a,  and  that  n = nQ,  T = T , at  s = 0,  this  model  (Adams 
and  Sturrock,  1975'  yields  the  following  expressions: 


T (s  ) = (l/  ''  - bFs 


n (s ) = 


(T0/T(S>)  exp  j-  f [(if  - bFs)5/?  - t5«  j , 


where  a 10  b ss  10  and  F (erg  cm  s ^ ) is  the  downward 

heat  flux. 

The  electrical  conductivity,  in  modified  gaussian  units,  may  be 
derived  from  the  expression  given  by  Spitzer  (1962): 


O = g T' 


where  g as  10  Hence  we  find,  from  (3.11),  that  5 is  related  to 


1? 


s by 


5 = 


I 

4 bgF 


To  • (To/2  - bFs)'/7 


(4.13 


Our  model  is  then  completely  specified  by  the  choice  of  the  coronal 
temperature,  the  coronal  density,  the  coronal  heat  flux,  y,  the  energy 
corresponding  to  and  the  injected  energetic  electron  flux.  For 
the  coronal  parameters,  we  have  adopted  values  typical  of  the  corona 
above  an  active  region  (Noyes,  1971): 

T =“  3 X 106  K 

Q -7 

n =-  10'  cm  0 

F ^ 5 X 10^  erg  cm  ^ sec  ^ 


We  choose  ^ to  correspond  to  25  kev  and  we  choose  y = 2.5,  corresponding 
to  an  x-ray  spectral  index  of  ~ 3 _ 3>5  (Brown,  1975 ) > which  is  a 
reasonable  value  for  the  spectral  index  at  peak  emission  (Datlowe,  19 75 ) • 
The  fraction  of  the  beam  energy  deposited,  and  the  total  energy 

& U 

deposited  by  joule  heating  between  T = 3 X 10  K and  T = 3 X 10  K as 

a function  of  the  energetic  electron  flux,  are  displayed  in  Figure  1. 

19.5  0 

For  a flare  area  of  10  r , the  energetic  electron  flux  inferred 

17  -2  - 1 

from  a large  impulsive  x-ray  burst  corresponds  to  ~ 10  cm  sec 
(Hoyng  et  al.,  1975)*  Figure  2 illustrates  the  energy  deposition 
rate  due  to  joule  heating  as  a function  of  temperature  of  the  atmosphere 
for  this  injected  energetic  electron  flux.  The  ordinate  of  Figure  2 
is  the  time  required  to  raise  the  ambient  plasma  temperature  by 
10  7 K. 
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V.  Discussion 


"1 


For  the  adopted  model.  Figure  2 shows  that  the  steady-state  assumption 
is  not  justified  for  the  energetic  electron  flux  required  to  produce  the 
observed  x-ray  flux  by  thick  target  bremsstrahlung.  The  atmosphere, 
especially  the  dense  low-temperature  region,  will  be  heated  rapidly  to 

j 

a high  temperature.  Once  the  plasma  is  raised  to  a few  times  10  K, 
ohmic  losses  will  become  unimportant  since  the  joule  heating  rate  is 
inversely  proportional  to  the  electrical  conductivity.  Joule  heating 
is  a more  effective  mechanism  for  heating  the  plasma  than  collisional 
losses  from  the  energetic  electrons,  because  the  ohmic  losses  are 
caused  by  thermal  electrons  in  the  reverse  current  which  have  much 
shorter  mean  free  paths  than  the  energetic  electrons.  It  is  likely 
that  a larger  fraction  of  the  beam  energy  is  deposited  than  is  indicated 

in  Figure  1 for  two  reasons.  First,  the  constant  heat  flux  model  under- 

estimates the  extent  of  the  low  temperature  region  where  the  ohmic  losses 
are  largest.  Second,  if  the  drift  velocity  associated  with  the  reverse 
current  becomes  too  large,  instabilities  can  develop  and  reduce  the 
electrical  conductivity  of  the  plasma.  If  the  conduct:*''  ity  of  the 
plasma  is  low  enough,  the  electric  field  that  develops  may  be  large 

enough  to  effectively  prevent  the  stream  from  reaching  the  chromosphere. 

The  time  scale  for  heating  the  dense  low- temperature  plasma  is 
comparable  with  the  time  scales  characteristic  of  impulsive  x-ray  bursts 
fHoyng  et  al.,  1975 )•  It  is  possible  that  thermal  bremsstrahlung  from 
the  rapidly  heated  plasma  can  account  for  a significant  portion  of  the 
observed  impulsive  x-ray  flux.  Hence  this  mechanism  offers  an  explanation 
of  the  fact  that  some  flares  first  produce  high-energy  x-ray  emission 


lU 


near  the  top  of  a loop  rather  than  at  the  foot-points  of  the  loop 
(Brueckner,  1976).  Another  important  consequence  of  this  process  is 
that,  if  thermal  emission  can  account  for  a substantial  fraction  or 
the  impulsive  flux  up  to  ~ 50  kev,  then  the  number  of  electrons 
required  to  produce  the  non-thermal  x-ray  flux  is  greatly  reduced 
[Brown,  1975). 

The  rapid  heating  of  the  plasma  will  also  cause  a large  pressure 
imbalance.  The  time  t (sec.)  for  the  plasma  to  respond  to  pressure 
changes  by  bulk  motions  can  be  estimated  from 


s 


where  L is  a characteristic  length  and  Vg  is  the  ion  thermal  velocity. 

I 2 

Even  for  a temperature  as  large  as  10  K,  this  time  is  long  (10  sec 

9*  7 

compared  to  the  heating  time  for  a length  scale  of  10^ * cm,  so  that 

the  plasma  density  will  not  change  appreciably  during  the  heating.  In 

4 

the  constant  heat  flux  model  we  have  adopted,  the  density  at  3 X 10  K 

is  approximately  100  times  the  density  at  3 X 10^  K.  If  both  these 

7 

regions  are  heated  to  a temperature  of  a few  times  10  K,  the  resulting 
pressure  imbalance  must  be  expected  to  produce  explosive  evaporation 
of  high  density  plasma  into  the  corona. 

A more  realistic  time-dependent  model  is  necessary  to  account 
for  the  complicated  phenomena  observed  in  solar  flares.  It  is  also 
possible  that  collective  effects  are  important  and  that  the  assumption 
of  a Maxwellian  velocity  distribution  for  the  ambient  plasma  must  be 
abandoned.  However,  our  analysis  of  reverse  current  in  terms  of  an 
idealized  model  strongly  suggests  that  joule  heating  is  likely  to  play 
an  important  role  in  the  heating  of  solar  plasma  during  a flare. 
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Figure  1.  The  fraction  of  the  beam  energy  deposited  (solid  curve  and 


the  total  energy  deposited  (broken  curve)  by  joule  heating 
as  a function  of  the  energetic  particle  number  flux. 

Figure  2.  Time  t required  for  ohmic  losses  to  heat  the  plasma  by 
10^  K as  a function  of  the  ambient  plasma  temperature. 
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